Proximity to a Nearly Superconducting Quantum Critical Liquid 
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The coupling between superconductors and a quantum critical liquid that is nearly superconducting 
provides natural interpretation for the Josephson effect over unexpectedly long junctions, and the 
remarkable stripe-spacing dependence of the critical temperature in LSCO and YBCO supercon- 
ductors. 
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In the resonating-valence-bond (RVB) theory, R] the 
superconductivity in underdoped cuprates is therniody- 
namically close to a quantum spin liquid. Although at 
zero doping there is a Neel ordered state, it is believed 
that after sufficient doping the mobilized holes disorder 
the Neel state to a quantum spin liquid. This picture is in 
general agreement with the experimental data measured 
in underdoped samples. For example, a quasi-particle 
pseudo-gap forms much above the superconducting tran- 
sition temperature Tc, and both the pseudo-gap and the 
energy gap below T^ exhibit the dj.2_y2-wave symmetry 
in the momentum space. \A Indeed, in a mean field solu- 
tion to the RVB theory, ^ the pseudo-gap follows from 
preformed spin pairs, and the condensation of these pairs, 
which is determined by the hole density x, defines the su- 
perconducting transition. The mean field solution cap- 
tures the major feature of the underdoped superconduc- 
tors, but it was accepted with skepticism because of the 
neglected gauge fluctuations. Q However, it was recently 
realized that the gauge fluctuations can be integrated out 
rigorously in the long wave length and low energy limit 
and a renormalized theory can be obtained, [pj In this 
respect, the notion of a quantum spin liquid is pertinent 
in the underdoped regime. 

In the cuprates, superconductivity is developed beyond 
a critical doping Xc{> xn), and the transition tempera- 
ture Tc increases with x at x > Xc (but x < Xp with Xp 
the optimum doping). Since the transition involves con- 
densation of spin pairs, it is reasonable to believe that 
the spin liquid at Xc is a quantum critical liquid (QCL) 
that is nearly superconducting, [g] The transition to a 
hole-uniform superconductor is, however, hidden by the 
phase separation instability at higher x, where the sys- 
tem is likely in a mixed phase composing of hole-rich 
stripes and the hole-poor critical liquid. |0] Below we ar- 
gue that proximity to the postulated QCL is consistent 
with two seemingly unrelated sets of data: the Joseph- 
son effect in unusually lengthy junctions made of under- 
doped cuprates at x slightly lower than Xc, ||] and the 
stripe-spacing dependence of the transition temperature 
of underdoped superconductors. ||-|lj] 



We first consider the general behavior of the QCL. It is 
known that the Coulomb interaction is strong and long- 
ranged in cuprates, which frustrates the phase separation 
of holes to form stripes. [Of Therefore, we believe that 
the dynamical exponent of the QCL should be z == 1. 
p3[ An immediate consequence of this exponent is that 
the inverse of a length has the dimension of energy. On 
the other hand, since the spin pairs are well formed, the 
amplitude fluctuations of the pairs are massive and can 
be integrated out in principle. The important dynamics 
comes from the phase 9 degrees of freedom of the pairs. 
Because of the isotropy in space and (imaginary) time, a 
generalized Josephson scaling relation for the phase stiff- 
ness can be obtained as follows. The Green's function for 
6 is, in general, 



Giq) = q 
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YiqO, 



(1) 



where q — (k, w) is a three-vector (k is the spatial 
wave vector and u> is the Matsubara frequency) , 77 is the 
anomalous dimension, ^ is the correlation length, and fi- 
nally Y{() is a scaling function of (. The phase stiffness p 
(including the superfluid density and the compressibility) 
is determined by [n4| 



lim G{q) = ¥/pq^ 



(2) 



where b = |(v^e*^)| ~ C '^^'^ is the order parameter with 
f3 and v being the scaling exponent for the order param- 
eter and correlation length, respectively. Here x is the 
hole density in the QCL. Thus 
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where in the last step we have used the generalized hyper- 
scaling relation 2/3 = v{d+z — 2 + r]) with d = 2 being the 
spatial dimension in our case. Thus the desired Joseph- 
son scaling relation reads p ^ ^^^ ■ This result could also 
be obtained by a simple power counting as follows. Near 
criticality, one expects that 



p{d,ef 



1 



(4) 



where the integration is over a volume of ^''+^ due to the 
isotropy in space and time. The ansatz 9^0 ^ 1/^ gives 
immediately the Josephson scaling. The above analysis 
implies that both the superfluid density and the com- 
pressibility scale as energy and/or inverse length. This 
statement is important to discuss the finite-size scaling 
in the following. With the general behavior of the QCL 
in mind, let us now discuss the relevance of the QCL in 
the cuprates. 

(I) Photodoping experiment: In a very recent and ele- 
gant experiment, H superconducting wires are generated 
by photo-doping a film with x slightly lower than Xc- The 
unluminated region serves as a junction. The critical cur- 
rent Ic is measured as a function of the junction length 
s (in the direction of current flow). It turns out that a 
relation 



l/s 



(5) 



holds from s = 45nm up to s = lOOnm, a length much 
larger than the coherence length S,s is^ ^ — bnm) in the 
superconducting leads. This phenomenon is difficult to 
understood if the junction is a normal insulator or metal, 
as one would expect IcS oc e~*' ' where I = ^s for a normal 
insulator and I = vp/T for a ballistic metal. According 
to the experimental setup, I <C lOOnm for both cases. 
Igl Therefore, the junction may be in a special state of 
matter, which is arguably a QCL. We note that the ge- 
ometry of the experimental setup is two one-dimensional 
(ID) superconducting wires immersed in a QCL. This is 
an important aspect of the experiment. The Cooper pairs 
can tunnel through an extended area connecting the two 
wires. Because the junction length is the only relevant 
length scale, the Josephson energy scales as 



Ej = F{s/i)/s 



(6) 



with a scaling function F[s/^) if the junction is not ex- 
actly at criticality. It is expected that F{s/£) remains 
a constant at s ^ ^. The above result can also be un- 
derstood as follows. Since s is the only length scale, the 
superfluid density ps in the junction sacles as l/s and the 
transverse extension of the supercurrent flow is Xg '^ s, 
up to corrections from a scaling function. The Josephson 
energy scale is thus given by 

Ej^ f dxdyp,{Vef^F{s/^)/s, (7) 

JsX, 

where the integration is over the area sXs where super- 
current flows. We have used V6 ^ l/s for a rough es- 
timate, and have included correction from the finiteness 
of s/S, by the scaling function. Since Ej determines the 
critical current Ic, in suitable units we have 



hs = F(s/C). 



(8) 



At s < ^ we expect that F is roughly a constant, so that 
Eq.(H) is recovered. The plot of /cS is nothing but a scal- 
ing function. From the data IcS drops to zero quickly 



at s > lOOnm. g This suggests that ^ ~ lOOnm. The 
decrease of IcS may be related to the thermal activation 
of vortices in the junction. More efforts should be de- 
voted to explain the sharp drop of the scaling function, 
which is beyond the scope of this scaling argument. On 
the other hand, because z = 1, we expect that ^ ^ 1/T 
in the critical regime, so that the scaling relation can be 
rewritten as 



IcS^F{Ts). 



(9) 



An important prediction follows from this function: In- 
stead of varying s, one could vary the temperature T to 
find the scaling behavior, provided that the superconduc- 
tivity within the wires is hardly changed. This situation 
is perhaps not feasible by photodoping, but is plausible 
by other techniques. 

To appreciate the importance of the special geome- 
try of the junction we addressed so far, consider another 
situation. Suppose the superconducting leads are wide 
enough so that effectively the QCL only lives within the 
junction. Neglecting boundary corrections the Josephson 
energy should be extensive in the junction width L. The 
only other length scale is s. By dimension analysis, we 
obtain 



dy dxp,{Vef^{L/s)F{s/ts/L)/s. (10) 



Ej 



In the last expression, the first factor L/s is extensive 
in L but is dimensionless. The remaining part is the 
characteristic energy. Here F{s/£,, s/L) is another scaling 
function that is finite at vanishing arguments. In this 
geometry, one would expect that 



IcS={L/s)F{s/ts/L). 



(11) 



At s <C ^ and s ^ L, one would expect that Ic oc L/s'^, 
which should be contrasted with Eq.(P). This scaling 
form leads us to propose the experimental test of the 
idea of QCL in the junction. If one can create a wide 
junction by photodoping the wide stripes of material then 
the critical current should follow Eq. (|ll|) . The data of 
Decca et al ^ begin with s > L, and are better described 
by the first geometry. (In fact, in their case, the width L 
can not be well defined as the wires become rounded at 
the tips). 

The QCL should also manifest itself in the magnetic 
field dependence of the Josephson critical current. Since 
the proliferation of vortices, each carrying a flux quan- 
tum, would wash out superconductivity, we expect that 
the magnetic field H defines another length scale through 
l/yH. Thus in the presence of a weak magnetic field, we 
expect that in a given geometry of the junction the field 
dependence of the critical current should be qualitatively 
given by 



Ic{H)/Ic{0) = 1 - ai^VH) 



where n is a universal positive exponent and a a non- 
universal positive constant. In the critical regime where 
^ > 1/T the above expression should be replaced by 
/c(i?)//c(0) = 1 - a{^/H/T^)"-. From the experiment 
of Decca et al, it turns out that n — 4. In order to 
see whether their results support the QCL scenario, it 
is desirable to vary both H and T and to see whether 
Ic{H)/Ic{0) is a scaling function of H/T^ (in the critical 
regime) . 

(II) Stripes in underdoped cuprates: Let us now dis- 
cuss the relevance of QCL in underdoped cuprates with 
hole-rich stripes. It was found that in LSCO materials 
there exists a striking linear relation between the super- 
conducting critical temperature T^ the incommensurate 
width S near the peak at (tt, tt) in the neutron scattering 
spectrum. |9|] Recently, a similar relation was discovered 
from the inelastic neutron scattering dada in YBCO sam- 
ples. ||l^ More interestingly, it was further argued that 
Tf oa S translates to 



Tc = hv* /s. 



(12) 



where h = 2TTh is the Planck constant, s the spacing be- 
tween stripes, and v* a material-dependent constant with 
the dimension of velocity. [^ Two aspects of Eq. (|l^) are 
quite unusual. First, for a given material, v* should be 
constant and small {v* ^ vp with vp being the Fermi 
velocity). And second, the robustness of the scahng up 
to the optimum doping, where the stripes, if they ex- 
isted at all, would be densely spaced (with s ^ Inm). It 
was conjectured that v* may be related to an unknown 
collective-mode of the stripes |lQ] and/or the motion of 
heavy quasiparticles near the flat band around (tt, 0). |l^] 
Given the experimental results, this conjecture is an im- 
portant clue to the underlying mechanism. Here we ar- 
gue that Eq. (O) follows from dynamic stripes coupled by 
QCL. 

We note that the stripes are metalic. They behaves 
as anti-phase domain walls for the remaining anti-ferro- 
magnetic background. |1^] (Recent mean-field calcula- 
tions for the t-J model indicated that in-phase stripes 
are slightly more favorable energetically. However, it 
was believed that quantum fluctuations might stablize 
anti-phase stripes. ||l^) The holes in the stripes are not 
strickly localized in the transverse direction. Rather, be- 
cause of quantum fluctuations, the stripes are dynamic 
unless pinned by inhomogeneities. Such transverse fluc- 
tuations were believed to develop strong superconduct- 
ing correlations along the longitudinal direction of the 
stripes. Q Also because of the transverse fluctuations 
of the stripes, the intervening region between the stripes 
may be driven to a quantum critical liquid, which is al- 
most superconducting. Combining these considerations, 
we model the whole system as a superposition of super- 
conducting stripes and quantum critical liquids plus cou- 
pling between them. pSl] This phenomenological model 



should be understood as an effective theory in the low 
energy limit. It does not answer why there are stripes 
and a QCL, which are the issue of a higher energy scale. 
Without loss of generality, let us assume that an isolated 
stripe a can be described by an action 



Dq 



1 



dydr [u{dy9„y + vidrOay 



where u and v are the bare stiffness components with 
respect to space y and time r, and are assumed to be 
independent of the the stripe spacing s. (Because of 
transverse fluctuations of the stripes, the spacing is only 
defined in an average sense.) As we argued, the stripes 
are coupled by the QCL. We imagine to find the effec- 
tive coupling between the stripes mediated by the QCL. 
This is virtually a real-space renormalization. In the long 
wave length and low energy limit, we need only to retain 
lowest order couplings. Clearly, there may be corrections 
to the bare stiffness components, which we believe is rela- 
tively small and is neglected in the following. Among the 
other effects, the Josephson-like coupling between near- 
est stripes a and (3, 



]^j dydTj{e^-epf 



is the most relevant. The coefficient j can be obtained by 
the same scaling arguments as that for the wide-junction 
geometry described above, 

where L now denotes the characteristic length of the 
stripes. Again in the long wave length and low energy 
limit, we can approximate {9a — OpY /s^ as {dxO)^ , so 
that the total action for the stripes can be put back into 
a continuum form in the r = (x, y, t) space as, 

s = ^Jd'rp^{d,ey, 

where /i denotes the three components, px = 
F{s/S,,s/ L)/ s, Py — u/s and p^- — v/s. This is an 
anisotropic 3D XY-model. To leading order all of the 
stiffness components in the new action are inversely pro- 
portional to s. If the system is macroscopically isotropic 
due to stripes in both directions, the global stiffness p is 
given by the geometrical mean of the three components, 
El so that 



p = H{s/i,s/L)/s. 



(13) 



Here H = F^/^^ is a scaling function given by F. The 
3D XY-model develops long range ordering as long as 
p > in thermodynamic limit (i.e., in an infinite plane 
and at zero temperature). [H This is again a quantum 
critical point with z = 1 due to isotropy in space and 
time. Since p has the dimension of energy, it naturally 



defines the scale of the superconducting gap. Thus at 
T > p the system is in a quantum critical state, while 
at T < p it is a renormalized superconductor. As a re- 
sult, the transiton temperature is given by Tc = p, which 
was also indicated earlier in some stripe scenarios of high 
temperature superconductivity. [Q In combination with 
the behavior of p we obtain 



T,s^H{s/C,s/L). 



(14) 



For long stripes, we can simplify the expression by the 
approximation TcS = H{s/^,0). This explains why Tc 
should be determined by s. The scaling function takes 
into account the correction due to the finiteness of s/£,. 
The above results follow essentially from a finite-size scal- 
ing with respect to a QCL. It is not clear a priori that 
such a scaling should work at all at the more or less mi- 
croscopic stripe spacing scales. Furthermore, the above 
relation does not guarantee that TcS will be a constant if 
H varies significantly at moderate arguments. However, 
^ ^ 1/Tc dX T = Tc in the critical regime of the QCL. 
Consequently, 



TcS = H{TcS,Qi). 



(15) 



This is a remarkable result, indicating that the plot of 
H(TcS, 0) against TcS only accesses a unique point of the 
scaling function at TcS = hv* , with hv* determined by 
the equation 



hv* = H{hv*,0). 



(16) 



We believe that this explains the robustness of Eq.(12) 



even though the scaling function H{(, 0) may vary with C 
due to the finite size effect. Clearly the notion of a QCL 
plays an important role in our arguments. 

In summary, proximity to a quantum critical liquid 
that is nearly superconducting provides appealing inter- 
pretation of the abnormal Joscphson effect and the stripe 
spacing dependence of Tc in underdoped cuprates. We 
also predicted the Joscphson effect in a geometry differ- 
ent from that in Ref. jg]. The field dependence of the 
Joscphson current is also discussed. Naturally one would 
expect other properties of superconductors in proximity 
to a quantum critical liquid. These may include scaling 
behavior of the various stiffness as a function of temper- 
ature, probing frequency and length scale. 
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